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Abstract 

The influence of vibrational motion on electron conduction through single molecules 
bound to metal electrodes is investigated employing first-principles electronic-structure 
calculations and projection-operator Green's function methods. Considering molec- 
ular junctions where a central phenyl ring is coupled via (alkane)thiol-bridges to gold 
electrodes, it is shown that - depending on the distance between the electronic tt- 
system and the metal — electronic-vibrational coupling may result in pronounced 
vibrational substructures in the transmittance, a significantly reduced current as 
well as a quenching of negative differential resistance effects. 



1 Introduction 



Recent advances in experimental studies of single molecule conduction [1,2,3,4,5,6,7] 
have stimulated great interest in the basic mechanisms which govern electron 
transport through nanoscale molecular junctions [8,9]. An interesting aspect 
that distinguishes molecular conductors from mesoscopic devices is the pos- 
sible influence of the nuclear degrees of freedom of the molecular bridge on 
electron transport. The current-induced excitation of the vibrations of the 
molecule may result in heating of the molecular bridge and possibly breakage 
of the junction. Conformational changes of the geometry of the conducting 
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molecule are possible mechanisms for switching behavior and negative dif- 
ferential resistance. Furthermore, the observation of vibrational structures in 
conduction measurements allows the unambiguous identification of the molec- 
ular character of the current. 



Vibrational structures in molecular conductance were observed, for example, 
in electron transport experiments on H 2 between platinum electrodes [3], as 
well as Cqo molecules between gold electrodes [2]. While in the latter two 
experiments the observed structures were attributed to the center of mass 
motion of the molecular bridge, other experiments on Ceo, C70 [7], and copper 
phtalocyanin [6] on an aluminum oxide film showed structures which were 
related to the internal vibrational modes of the molecular bridge. Moreover, 
vibrational signatures of molecular bridges have also been observed in inelastic 
electron tunneling spectroscopy (IETS) [10,11]. 



This tremendous experimental progress has inspired great interest in the 
theoretical modelling and simulation of vibrationally-coupled electron trans- 
port in molecular junctions. In the low-voltage, off-resonant transport regime 
combinations of electronic-structure calculations and nonequilibrium Green's 
function theory (NEGF), employing the self-consistent Born approximation 
(SCBA), have been used to investigate vibrational signatures in IETS [12,13,14]. 
The majority of the studies in the resonant tunneling regime (for higher volt- 
ages) have been based on generic tight-binding models, using the NEGF-SCBA 
approach [15], or kinetic rate equations to calculate the current [16,17,18,19,20,21]. 
These model studies have demonstrated that the vibrational motion of the 
molecular bridge may affect the current-voltage characteristics significantly. 



Most of the studies of vibrationally-coupled electron transport reported so 
far invoke approximations which restrict their applicability either to small 
electronic-vibrational coupling, small molecule-lead coupling or separability 
of the vibrational modes. To circumvent this limitation, we recently pro- 
posed an approach [22,23,24], which is based on scattering theory and the 
projection-operator formalism of resonant electron- molecule scattering [25]. 
Within the single-electron description of electron conduction, this approach 
is valid for strong electronic-vibrational and molecule-lead coupling and thus 
allows the study of electron transport in the resonant regime. In the work re- 
ported here, the approach is combined with first-principles electronic-structure 
calculations to characterize the molecular junction. It is applied to investigate 
electron transport in two systems: p-benzene-dithiolate (BDT) and p-benzene- 
di(ethanethiolate) (BDET) covalently bound to gold electrodes (Fig. 1). 
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Fig. 1. The two molecular junctions investigated: BDT (top) and BDET (bottom) 
between clusters of 30 gold atoms. 

2 Theory 

To study the influence of vibronic effects on single molecule conductance in 
BDT and BDET, we use an ab-initio based model described by the Hamilto- 
nian 

H= y £\<f> j )E j (<f> j \+ Yl \<l>k)E k (<f> k \+V + H n , (la) 

j&M keL,R 

v=Y E (l^VjM + mvMi). (ib) 

jSM k£L,R 

Here, \4>j), j £ M and \4> k ), k G L,R denote orthogonal electronic states 
(molecular orbitals) representing the molecular bridge and the left and right 
leads, respectively, Ej and E k are the corresponding energies, and V charac- 
terizes the coupling strength between molecule and leads. The nuclear degrees 
of freedom as well as the electron- nuclear coupling are described by H n . Be- 
cause the systems studied in this work do not exhibit large amplitude motion, 
we can employ the harmonic approximation and a low order expansion of 
the molecular orbital energies around the equilibrium geometry of the neutral 
molecule. Thus, the nuclear part of the Hamiltonian reads 



H n = H n0 + H ne = ^uiajai + \<t>j) K \ f zMjV ( 2 ) 
i jeM,l 

Here, aj and ai are creation and annihilation operators for the Z-th normal 
mode with dimensionless coordinate qi and frequency ui, and /Cj denotes the 
corresponding vibronic coupling constant in state \<j)j). 
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For both systems investigated, the parameters of the Hamiltonian (1) were 
determined employing electronic-structure calculations [26]. A detailed de- 
scription of the strategy used will be given in a future publication. Briefly, 
the overall system was first separated into molecule and leads. The molecule 
and those parts of the leads which are close to the molecule (in the following 
referred to as the contacts) where treated explicitly in the quantum chemical 
calculations. A finite cluster of 30 gold atoms on each side of the molecule 
with a tip-like geometry was employed to model the contacts (cf. Fig. 1). A 
partial geometry optimization of the molecule and the first layer of the gold 
clusters was employed to determine a realistic molecule-lead binding geome- 
try. Thereby, the sulfur atoms are bound covalently to two gold atoms, which 
is the preferred bond formation if no symmetry constraints are applied [4,27]. 
The effect of infinite leads can in principle be described employing the surface 
Green's function for the contacts [28]. In the present study, we have added 
constant selfenergies to the atomic orbital energies of the outer gold atoms 
of the contacts. After a Lowdin orthogonalization of the Kohn-Sham orbitals, 
the overall system was partitioned into molecule and contacts using projection 
operator techniques [29]. A subsequent partial diagonalization of the Hamil- 
tonian within the three subspaces (left contact, molecule, right contact) gives 
the electronic states \4>j) and \<f>k), representing the molecular bridge and the 
leads, the energies Ej, as well as the molecule- lead coupling matrix elements 
Vjk- The nuclear degrees of freedom of the molecular bridge were character- 
ized based on a normal mode analysis of an extended molecule that includes 
seven gold atoms on each side of the molecule. The electronic-nuclear coupling 
constants were obtained from the gradients of the electronic energies Ej 
with respect to the normal coordinates 

To investigate vibrationally inelastic electron transport through molecular 
junctions we shall consider the transmission probability of a single electron 
through the junction and the current- voltage characteristics. The inelastic 
transmission probability of a single electron from the left to the right lead as 
a function of initial and final electron energy is given by the expression 

T R ^ L (E h E f ) = 4tt 2 ]T ]T £ P Vi S(E f + E Vf -Ei- E Vi ) 

v i. v / ki£L k f £R 

x6(Ei- E ki )5(E f - E^KvfKkflVGiEjVlh)^)] 2 .(3) 

Here, the 5-function accounts for energy conservation Ef + E Vf = E^ + E Vt 
with E v . and E Vf being the energy of the initial and final vibrational states 
|vj), \vf), respectively, and P Vi = (vj|p |vj) denotes the population probability 
of the initial vibrational state po — e~ Hn0 ^ kBl " > / Z . In the systems considered 
in this work, the electron couples primarily to modes with relatively high 
frequencies. As a result, thermal effects are not expected to be of relevance 
and the initial vibrational state is assumed to be in the ground state, i.e. 
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Po = |0)(0|. The total transmission probability, T R ^i(Ei), is obtained by 
integrating T R ^ L (Ei, Ef) over the final energy of the electron. 

To calculate the transmission probability, the Green's function G(E) = {E + — 
H)^ 1 is projected onto the molecular space using the projection operators 

P = EjeM\<f>j)(<f>j\, Ql = EfeeL \(f>k)(<f>k\, and Q R = E fcG R |0fc) (0fc|, resulting 
in the expression 



T R ^ L (Ei, E f ) = £ tr {8(E f + E Vf - E t - H n0 )p (4) 
v / 

x T L (E t ) Gl(Ei) \v f )(v f \r R (E f ) G M (E t )}, 
with the Green's function projected on the molecular bridge 

G M (E) = PG(E)P = (E + - PHP - E L (E) - ^r(E))- 1 (5a) 
E L (E) = PVQ L (E + - Q L HQ L y x Q L VP = - % -Y L {E) + A L (E). (5b) 

Here, Hl(E) denotes the self energy due to coupling to the left lead and similar 
for the right lead. 

It is noted that in contrast to purely electronic transport calculations as well 
as applications of the non-equilibrium Green's function formalism to vibronic 
transport [12,30], the Green's function G M and the self energies S L , T, R are 
operators both with respect to the electronic and nuclear degrees of freedom. 
Thus, no separability of the nuclear degrees of freedom is assumed and the 
Green's function has to be evaluated in the combined electronic-vibrational 
Hilbert space. To reduce the computational effort, in the calculations presented 
below, the four vibrational modes with the strongest vibronic coupling (as de- 
termined by the ratio of the electronic-vibrational coupling and the electronic 
coupling) were explicitly taken into account. Those are the C-C-C bending, 
the ring breathing (only for BDT), the C-C-H bending, and a C-C stretching 
mode, depicted in Figs. 2, 3. The corresponding vibronic coupling constants in 
the most important molecular orbitals are given in Tables 1, 2. Furthermore, 
the number of states \<f>j) on the molecular bridge, which were explicitly in- 
cluded in the calculation of the transmission, was reduced by including only 
those with energies in the vicinity of the Fermi energy (7 for BDT and 12 for 
BDET). 

Based on the transmission probability (4), the current through the molecular 
junction is calculated using the generalized Landauer formula [31] 
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Fig. 2. Normal modes of BDT included in the calculation. 




Fig. 3. Normal modes of BDET included in the calculation. 

Table 1 

Frequencies of the four most important vibrational modes of BDT as well as gra- 
dients of these modes for the three orbitals (A,B,C) contributing to conduction. 





uj (cm 1 ) 


k a (10" 1 eV) 


K B (10" 1 eV) 


k c (10- 1 eV) 


(a) 


349.26 


0.61 


0.08 


0.62 


(b) 


1092.02 


0.72 


0.80 


0.33 


(c) 


1198.14 


0.90 


0.13 


0.16 


(d) 


1627.28 


1.52 


0.93 


0.36 



Table 2 

Frequencies of the four most important vibrational modes of BDET as well as 
gradients of these modes for the two orbitals (D,E) contributing to conduction. 





uj (cm 1 ) 


k d (10" 1 eV) 


k e (10" 1 eV) 


(a) 


544.48 


0.76 


0.22 


(b) 


1197.11 


0.51 


0.69 


(c) 


1229.49 


1.10 


0.47 


(d) 


1671.56 


1.36 


1.62 



6 



-T^niE^E^UE^l- f L {E f )]}, (6) 

where fi(E), fn{E) denote the Fermi distribution for the left and the right 
lead, respectively. While the formulas (3), (4) for the single-electron transmis- 
sion probability involve no approximation, the expression (6) for the current 
is valid if many-electron processes are negligible. Furthermore it is implicitly 
assumed, that the nuclear degrees of freedom of the molecular bridge relax to 
the vibrational equilibrium state after transmission of an electron. 

In principle, the basis states \4>j), \(j>k), the electronic energies and the nuclear 
parameters depend on the bias voltage. For the studies in this work, we have 
used, for simplicity, parameters obtained from electronic-structure calculations 
at equilibrium and assumed that the bias voltage V enters the formulas only 
via the chemical potentials of the leads ^l/r = e f =t eV/2. Here €f denotes 
the Fermi energy at equilibrium, which was approximated as the average of 
the energies of the HOMO and LUMO orbitals of the overall system . The 
energies of the lead states for finite voltage are thus given by E k ±eV/2. Since 
we do not invoke the wide-band approximation, the Green's function Gm, the 
self energies Sl/_r as well as the width functions Yl/r also depend on the bias 
voltage. 



3 Results and Discussion 

We first consider electron transport through BDT. Fig. 4a shows the to- 
tal transmission probability for this system. In addition to the transmission 
based on a vibronic calculation, also the result of a purely electronic calcula- 
tion (where all electronic-nuclear coupling constants were set to zero) is 
shown. The transmission probability exhibits several pronounced peaks. The 
first three peaks below e/ at energies -2.61, -2.06, and -1.25 eV, respectively, 
are caused by the three orbitals depicted in Fig. 4a. These orbitals are related 
to the ei 9 -orbitals of benzene and p-orbitals at the sulfur atoms. While orbitals 
B and C have significant contributions from the bridging sulfur atoms, orbital 
A is completely localized on the phenyl ring. As a result, the molecule- lead 
coupling strength for the three orbitals is quite different: the width function, 
Yjj, at the respective peak energy varies from 0.01 eV (A), 0.07 eV (B) to 
0.25 eV (C), respectively. Consequently, the structures in the transmission 
probability caused by orbitals B and C are rather broad, whereas orbital A 
results in a narrow peak. Besides structures at higher energies, the transmis- 
sion probability also exhibits a broad peak close to e/. Although this peak 
is influenced by orbitals B and C, it cannot be directly related to any of the 
orbitals of the molecular bridge and is therefore attributed to a metal-induced 
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Fig. 4. (a): Electronic (dashed line) and vibronic (solid line) total transmission prob- 
ability through a BDT molecular junction at zero voltage as a function of the initial 
energy of the electron (relative to the Fermi energy) . The three localized molecular 
orbitals, denoted A, B, and C, dominate the transmittance at the indicated peaks, 
(b): Current- voltage characteristic of BDT. Shown are results of calculations with 
(solid line) and without (crosses) coupling to molecular vibrations. 

gap state [32,28]. As is known from scanning tunneling spectroscopy [33], the 
appearance of such states is expected for tip-like geometries of metal contacts 
as used in the present study. 

The comparison of the vibronic transmission probability (solid line in Fig. 
4a) with the results of the purely electronic calculation (dashed line) reveals 
that the electron-nuclear coupling in BDT has significant effects on the trans- 
mittance for narrow resonances. In particular, it results in a splitting of the 
narrow peak at -2.61 eV into a number of smaller structures. The effect of 
nuclear motion on the other peaks in the transmission probability is, on the 
other hand, rather small. This can be rationalized by considering the vibronic 
and electronic coupling constants in the corresponding states. The importance 
of vibronic effects caused by the nuclear mode qi in state \(f>j) is determined by 
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the ratio of the vibronic and electronic coupling constants, Kf' /Yjj. Although 
the electronic- vibrational coupling constants for some of the nuclear modes 
in states \<Pb) and \<pc) are relatively large (cf. Table 1), the lifetime of the 
electron on the molecular bridge (approximately given by h/Tjj) is short - 
thus resulting in a small effective coupling /Yjj. 

The simulated current-voltage characteristic of BDT, depicted in Fig. 4b, ex- 
hibits a nonlinear behavior: A small increase of the current at small voltage 
caused by the metal-induced gap state at e/ is followed by a strong increase at 
larger voltage resulting from the contributions of orbitals B and C. Because 
orbitals B and C are strongly coupled to the metal leads and thus have a small 
effective vibronic coupling, the influence of nuclear motion on the current in 
BDT is almost negligible. 

We next consider electron transport through BDET. In BDET, the (CH 2 ) 2 - 
spacer groups reduce the electronic coupling between the conjugated 7r-system 
of the phenyl ring and the metal leads, which has a number of interesting 
consequences for the transport mechanism. The transmission probability for 
BDET is depicted in Fig. 5a. In contrast to the result for BDT, the transmis- 
sion probability in the BDET system has only very little intensity at e/. Thus, 
metal-induced gap states are of minor importance in this system. This is a con- 
sequence of the larger size of BDET and the smaller electronic coupling [28]. 
The resonance peaks at energies -2.24 eV and -1.85 eV, which are closest to 
ej and thus determine the transport process, are caused by molecular orbitals 
D and E. While orbital D resembles an ei 9 -orbital of benzene, orbital E has 
additional contributions at the spacer groups and the sulfur atoms. As a con- 
sequence, state D has small electronic coupling to the leads (Y DD — 2.1 • 1CT 4 
eV) resulting in a very narrow peak in the transmission probability, whereas 
the significant coupling of orbital E to the leads (Tee = 8.9-10 -2 eV) results in 
a rather broad structure. Compared to the corresponding orbitals of the BDT 
system, the spacer group reduces the molecule-lead coupling by about an order 
of magnitude. As a result, the lifetime of the electron on the molecular bridge 
is significantly longer and thus effects due to nuclear motion become more 
important. The comparison between the results of vibronic (solid line in Fig. 
5a) and purely electronic (dashed line) calculations demonstrates that the 
electronic-vibrational coupling in BDET alters the transmission probability 
significantly. In particular, the narrow peak (D) splits into many vibrational 
subpeaks, but also the broader peak (E) acquires substantial vibrational sub- 
structure. The vibrational structures can be assigned to the individual nuclear 
modes of BDET [29]. 

The current-voltage characteristics of BDET is shown in Fig. 5b. The result 
obtained from a purely electronic calculation (dashed line) shows an increase 
of the current at about 3.5 Volt caused by state E, which is followed by a pro- 
nounced decrease of the current at 4 Volt. The weakly coupled state D results 
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Fig. 5. (a): Electronic (dashed line) and vibronic (solid line) total transmission 
probability through a BDET molecular junction at zero voltage as a function of 
the initial energy of the electron (relative to the Fermi energy). The two orbitals, 
denoted D and E, dominate the transmittance at the indicated peaks, (b): Cur- 
rent-voltage characteristic of BDET. Shown are results of calculations with (solid 
line) and without (dashed line) coupling to molecular vibrations. 

only in a small step-like increase of the current at about 4.5 Volt. A detailed 
analysis [29] shows that the negative-differential resistance (NDR) effect at 
4 Volt is caused by the voltage dependence of the self energies Sx, and 
the corresponding width functions T. It should be emphasized that this NDR 
effect can only be described if the energy dependence of the self-energies is 
taken into account and will be missed within the often used wide-band ap- 
proximation. Including the coupling to the nuclear degrees of freedom changes 
the current-voltage characteristics substantially. In particular, the electronic- 
vibrational coupling results in a quenching of the NDR effect, a significantly 
smaller current for voltages in the range 3.5 - 4.25 V, and noticeable vibra- 
tional structures in the current. 

Finally, it is noted that the overall current through BDET is about an order of 



10 



magnitude smaller than in BDT. This is a result of the (CH^-spacer groups, 
which reduce the effective coupling between phenyl-ring and gold leads. A 
similar reduction of the current was also found in experimental studies of 
electron transport through benzenedimethanethiol [5]. 



4 Conclusion 

The theoretical studies reported in this work demonstrate that nuclear motion 
can have a significant effect on electron conduction through single molecule 
junctions. The importance of vibronic effects thereby depends crucially on 
the relative ratio between electronic-vibrational coupling and molecule-lead 
interaction. While the former is determined by the character of the electronic 
states of the molecular bridge, the latter can be varied systematically, e.g., by 
introducing spacer groups between molecule and leads. To study these effects, 
we have considered in this work two systems, BDT and BDET, covalently 
bound to gold electrodes. In BDT, an example for a molecular junction with 
strong coupling to the leads, electronic-vibrational coupling results in vibra- 
tional sidebands in the transmission probability, but the effect is rather small 
and the influence on the current is almost negligible. In BDET, on the other 
hand, a system where spacer groups reduce the molecule-lead coupling by 
about an order of magnitude, vibronic effects change the current- volt age char- 
acteristics substantially. In particular, it was found to result in quenching of 
NDR-effects, a significantly reduced current over a range of voltages, as well 
as vibronic structures in the current- volt age characteristics. 



5 Acknowledgment 

This work has been supported by the Deutsche Forschungsgemeinschaft. 



References 

[1] M. Reed, C. Zhou, C. Muller, T. Burgin, J. Tour, Science 278 (1997) 252. 

[2] H. Park, J. Park, A. Lim, E. Anderson, A. Alivisatos, P. McEuen, Nature 
(London) 407 (2000) 57. 

[3] R. Smit, Y. Noat, C. Untiedt, N. Lang, M. van Hemert, J. van Ruitenbeek, 
Nature (London) 419 (2002) 906. 



11 



[4] H. B. Weber, J. Reichert, F. Weigend, R. Ochs, D. Beckmann, M. Mayor, 
R. Ahlrichs, H. v.Ldhneysen, Electronic transport through single conjugated 
molecules, Chem. Phys. 281 (2002) 113-125. 

[5] X. Xiao, B. Xu, N. Tao, Nano Lett. 4 (2004) 267. 

[6] X. Qiu, G. Nazin, W. Ho, Phys. Rev. Lett. 92 (2004) 206102. 

[7] N. Liu, N. Pradhan, W. Ho, J. Chem. Phys. 120 (2004) 11371. 

[8] P. Hanggi, M. Ratner, S. Yaliraki, (Eds.), Chem. Phys. special issue on: 
"Processes in molecular wires" 281. 

[9] A. Nitzan, M. Ratner, Science 300 (2003) 1384. 

[10] B. Stipe, M. Rezai, W. Ho, Science 280 (1998) 1732. 

[11] J. Kushmerick, J. Lazorcik, C. Patterson, R. Shashidhar, D. S. Seferos, G. C. 
Bazan, Nano Lett. 4 (2004) 639. 

[12] A. Pecchia, A. D. Carlo, Nano Lett. 4 (2004) 2109. 

[13] Y. Chen, M. Zwolak, M. D. Ventra, Inelastic current-voltage characteristics of 
atomic and molecular junctions, Nano Lett. (5) (2005) 813. 

[14] M. Paulsson, T. Frederiksen, M. Brandbyge, Phys. Rev. B 72 (2005) 201101. 

[15] M. Galperin, M. Ratner, A. Nitzan, Phys. Rev. B 73 (2006) 045314. 

[16] E. Emberly, G. Kirczenow, Phys. Rev. B 61 (2000) 5740. 

[17] D. Boese, H. Schoeller, Europhys. Lett. 54 (2001) 668. 

[18] V. May, Phys. Rev. B 66 (2002) 245411. 

[19] J. Lehmann, S. Kohler, V. May, P. Hanggi, J. Chem. Phys. 121 (2004) 2278. 

[20] J. Koch, F. von Oppen, Phys. Rev. Lett. 94 (2005) 206804. 

[21] K. Nowack, M. Wegewijs, New J. Phys. 7 (2005) 239. 

[22] M. Cizek, M. Thoss, W. Domcke, Phys. Rev. B 70 (2004) 125406. 

[23] M. Cizek, M. Thoss, W. Domcke, Czech. J. Phys. 55 (2005) 189. 

[24] A closely related approach has been proposed by Ness et al. [34]. 

[25] W. Domcke, Phys. Rep. 208 (1991) 97. 

[26] All quantum chemical calculations were performed with the TURBOMOLE 
package (V5-7) [35] using density fuctional theory (DFT) with the B3-LYP 
functional and the SV(P) basis set including ECP-60-MWB on the gold atoms. 

[27] H. Basch, M. A. Ratner, Binding at molecule/gold transport interfaces I. 
geometry and bonding, J. Chem. Phys. 119 (2003) 11926. 

[28] Y. Xue, M. Ratner, Phys. Rev. B 68 (2001) 115406. 



12 



[29] C. Benesch, M. Cizek, W. Domcke, M. Thoss, to be published. 

[30] D. A. Ryndyk, M. Hartung, G. Cuniberti, Phys. Rev. B 73 (2006) 045420. 

[31] A. Nitzan, Electron transmission through molecules and molecular interfaces, 
Annu. Rev. Phys. Chem. 52 (2001) 681-750. 

[32] R. Guitierrez, G. Fagas, K. Richter, F. Grossmann, R. Schmidt, Europhys. Lett. 
62 (2003) 90. 

[33] A. V. de Parga, O. Hernan, R. Miranda, A. L. Yeyati, N. Mingo, A. Martin- 
Rodero, F. Flores, Phys. Rev. Lett. 80 (1998) 357. 

[34] H. Ness, S. Shevlin, A. Fisher, Phys. Rev. B 63 (2001) 125422. 

[35] R. Ahlrichs, M. Bar, M. Haser, H. Horn, C. Kolmel, Chem. Phys. Lett. 162 
(1989) 165. 



13 



